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VERTICAL  RESOLUTION  OF  TEMPERATURE  PROFILES  FOR 
HIGH  RESOLUTION  INFRARED  RADIATION  SOUNDER  (HIRS) 


Y.  M.  Chen",  H.  M.  Woolf  and  W.  L.  Smith 
National  Environmental  Satellite  Service 
National  Oceanic  and  Atmospheric  Administration 
Washington,  D.  C. 


ABSTRACT.   By  using  the  method  of  Backus  and  Gilbert, 
the  intrinsic  vertical  resolution  of  temperature  pro- 
files obtained  from  the  High  Resolution  Infrared 
Radiation  Sounder   (HIRS)  on  the  Nimbus-F  satellite, 
which  contains  channels  in  both  4.3-ym  and  15-ym  CO2 
absorption  bands,  is  evaluated  and  compared  with  those 
of  channels  in  the  M-.3-ym  band  and  the  15-um  band  alone, 
It  is  found  that  the  combination  of  4.3-ym  and  15-ym 
bands  is  superior  to  the  other  two  cases  for  all  levels 
of  the  pressure. 


1.   INTRODUCTION 

One  of  the  basic  problems  in  atmospheric  temperature  profile  inversion 
techniques  is  the  intrinsic  vertical  resolution  of  the  measurement  system. 
Recently,  the  fundamental  technique  for  obtaining  a  direct  quantitative 
value  in  answer  to  the  problem  was  developed  by  Backus  and  Gilbert  (1968, 
1970)  for  the  analysis  of  the  solid  earth.   Conrath  (1971)  has  applied  it 
to  the  problem  of  the  vertical  sounding  of  the  atmosphere  by  means  of 
remote  radiation  measurements  (with  random  noise)  in  the  15-ym  CO2  absorp- 
tion band. 

The  purposes  of  this  paper  are  to  evaluate  the  intrinsic  vertical  resolu- 
tion of  temperature  profiles  obtained  from  the  High  Resolution  Infrared 
Sounder   (HIRS)     on  the  Nimbus-F  satellite  (which  contains  channels  in 
both  4.3-ym  and  15-ym  CO2  absorption  bands)  and  to  compare  it  with  those 
obtained  separately  in  the  4-,3-um  and  15-ym  bands.   First,  a  brief  review 
of  the  formulation  of  Backus  and  Gilbert  is  presented.   Then  for  various 
levels  of  atmospheric  pressure  p,  the  averaging  kernel  A(x,x'),  the  center 
of  A(x,x'),  and  the  resolving  length  of  A(x,xT)  are  computed  and  plotted 
as  functions  of  p  for  three  individual  cases,  the  4.3-ym  band,  the  15-ym 
band  and  the  4.3-ym  +  15-ym  bands  (HIRS).       Finally,  we  discuss  and 
compare  these  three  cases. 

2.   THE  METHOD  OF  BACKUS  AND  GILBERT 

In  the  remote  sensing  of  atmospheric  temperature  profiles,  measurements 
are  made  of  radiances  in  a  finite  number  of  spectral  intervals  within 
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atmospheric  absorption  bands.   For  a  nonscattering  atmosphere  in  local 
thermodynamic  equilibrium,  the  expression  for  the  radiances  has  the  form 

r  i  rx(Po)      r 

I(vi)  =  B  T(Po),vi  T(Po,vi)  -\  B   T(p),vi"|  d^^v^  dx(p)  ,    (l) 

J  Jx(o)     L       -1  dx^P) 

i  =  1j  2, . . . ,  m  , 

where  B[T(p),v]  is  the  Planck  function  for  temperature  T(p)  at  pressure 
p;  the  independent  variable  x  can  be  any  monotonic  function  of  p  (x  =  In  p 
being  used  here);  pQ  is  the  surface  pressure;  and  x(p9uj_)  is  the  trans- 
mittance  of  the  atmosphere  above  pressure  p  at  u:. 

The  inversion  problem  then  is  to  estimate  the  temperature  profile  T(p), 
given  the  atmospheric  transmittances  x(p,v. )  and  measurements  of  the 
radiances  (l(vi)},  i  =  1,  2,...,  m.   If  the  contribution  from  the  non- 
linearity  of  eq  (1)  is  small,  then  only  the  first  iteration  of  the  iterative 
algorithm  of  Chen  (1973)  or  the  procedure  of  Backus  and  Gilbert  will  be 
sufficient  for  our  purpose.   Hence,  the  contribution  from  nonlinearity  is 
assumed  to  be  small  throughout. 

Let 

T  (p)  =  T0(p)  +  6T(p)  ,  (2) 

where  T0(p)  is  the  reference  profile  and  is  taken  to  be  the  U.S.  Standard 
Atmosphere.   Then  from  Chen  (1973),  6T(p)  should  be  the  best  approximate 
solution  of  the  linearized  equation 


fxo 

I     kIx.v^   6T(x)   dx   =   6I(v1)  ,    i  =  1,  2,..., 


m.  (3) 


where      K  [x,v.  |  =   d  B  [t(x),  v-]  i    dx(x,v1- ) 

L   1-l         d   T        T     dx 

I  o 

$x 
fi|T0(x)>vil  dx(x,vn-)  dx  -  I(vi) 


(4) 


dx  (5) 


and  dB/dT  is  the  Frechet  derivative  of  b[t], 


The  weighted  average  of  61   at  x  which  gives  heavy  emphasis  to  points  close 
to  x  and  very  little  to  distant  points  is  defined  by 

Sxo 
A(x,x')   6T(x')   dx'  ,  (6) 

_00 

where  the  averaging  kernel  A  is  normalized  according  to 


\ 


xo 


A(x.x')   dx1   =  1  .  (7) 


<6T>   is  the  most  localized  weighted  average  if  and  only  if  the  averaging 

kernel  A  closely  resembles  the  Dirac  delta  function  5(x'-x).   The  spread  of 
h.   from  x  is  defined  by 

Q(  x,  A  )  =  aT\  °  J(x,x')  A2(x,x')  dx'  ,  (8) 


where  J(x,x?)  is  a  chosen  infinitely  differentiable  function  of  x'  such  that 
J(x,x)  =  0  and  increases  monotonically  as  x'  increases  or  decreases  away 
from  x  with  dimension  of  x'^,  and 


aJ 


i: 


2  (9) 


J(x,x')  Aj(x,x')  dx' 


with    A£(x,x')  =  l/l    ,   x  -  */2  <x'  <x  +  Z/2    , 


(10) 


0,   |x'  -  x|^A/2  , 

Note  that    Q(  x,  Ajj  )  =  I        (  the  spread  of  Aj  ).  (11) 

If  A  has  a  sharp  peak,  the  choice  of  J  is  rather  free.   However,  if  A  has  a 
blurred  peak,  the  selection  of  J  is  much  more  critical.   Throughout, 

J(x,x')  =  (  x  -  x'  )2  (12) 

is  used.   (otj  =  12.)   As  one  will  see  later,  this  choice  of  J  is  not  very 
good  for  4.3-ym  CO2  channels  in  the  upper  stratosphere  (   p  =  8  -  60  mb). 

Equation  (3)  shows  that  {6I(v-)}  ,  i  =  l,2,...,m,  are  m  bounded  linear 
functionals  of  6T(x).   Since  <6T>   and  {6I(vi)}J  i  =  l,2,...,m,  all  de- 
pend linearly  on  the  function  6T,  it  follows  that  <6T>X  must  depend  line- 
arly on  {Sl(v^)}.   Therefore,  there  should  exist  constants  {aj'(x)},  i  =  1,2, 
. . . ,m,  depending  on  the  fixed  point  x  such  that 

<6T>x  =  £    ai<x>  6i<vi>  •  (i3> 

1=1 

Hence  from  (3)  and  (6), 

A(x,x')   =   J    ai(x)   K|x',Vil  .  (14) 

It  follows  from  (8)  and  (14)  that  Q(x,  A)  is  a  positive-definite  quadratic 
function  of  {a^(x)}.   Since  the  determination  of  the  most  localized  weighted 


average  <5T>   is  equivalent  to  minimizing  the  spread  of  A  from  x,  subject 
to  the  constraint  (7),  by  using  the  method  of  Lagrange  multipliers,  the 
proper  set  of  {a^(x)}  satisfies  the  following  set  of  m+1  linear  algebraic 
equations : 


in 

5=1 


1 12  I    "    (x-xM2K|~xf  ,v  IkTx'.v.      dx'l    a.(x)   +  Al    ^fx'^     dx'=0, 


ki!-l*^] 


i   =   1,2,. . . ,m. 


(15) 


dx'J>     a.(x)   =   1    , 


where  A  is  a  Lagrange  multiplier. 

The  vertical  resolution  of  the  temperature  profile  at  the  level  x  obtain- 
able from  a  given  set  of  radiance  measurements  is  determined  by  the  close- 
ness of  A(x,x')  to  6(x'  -  x).   Visually,  the  plot  of  A  (x,xT)  vs.  x'  gives 
a  qualitative  estimate  of  the  vertical  resolution  at  x.   To  characterize 
the  behavior  of  A(x,x')  more  precisely,  the  "resolving  length"  of  A(x,x') 
is  introduced  and  defined  as  the  spread  about  its  center, 


w(x)  =  12 


\!  [c(x)-x] 


c(x)  -  x'l  2  A2(x,x')  dx'  , 


(16) 


where 


fXo      2  /fXo    O 

c(x)  =  I    x1  A^(x,x'>  dx'/  1     A2(x,x')  dx' 


(17) 


is  the  "center"  of  A(x,x'). 


3.   APPLICATION  TO,  HIRS 

The  method  of  Backus  and  Gilbert  has  been  applied  to  an  analysis  of  the 
vertical  resolution  of  temperature  profiles  inferred  from  radiance  measure- 
ments in  two  CO2  absorption  bands  observed  by  the  High  Resolution  Infrared 
Radiation  Sounder  (HIRS) ,    The  7  channels  in  the  15-um  C09  band  and  the 
5  channels  in  the  4.3-um  CO2  band  in  the  study  are  listed  in  table  1.   The 
dx^/dx(p)  vs.  p  curves  for  15-um  and  4.3-um  C0?  bands  are  shown  in  figures 
1  and  2  respectively.   Similarly,  the  dB^/dT |T_.dT£/dx(p)  radiative  transfer 
kernels  vs.  p  for  15-um  and  4.3-um  CO2  bands  are  shown  in  figure  3.   The 
averaging  kernel  A  for  each  of  various  pressure  levels,  ranging  from  1  mb 
to  1000  mb,  is  computed  for  3  separate  cases:   the  5  channels  in  the  4-,3-um 
band  alone,  the  7  channels  in  the  15-um  band  alone,  and  the  12  channels  in 
the  combination  of  4.3-um  and  15-um  bands.   The  corresponding  set  of  three 
averaging  kernels  at  various  levels  of  p  are  plotted  as  functions  of  atmos- 
pheric pressure  in  figures  M-  to  12.   These  curves  give  a  qualitative  esti- 
mate of  the  vertical  resolution  at  the  indicated  levels  and  offer  a  visual 
comparison  for  the  resolving  length  and  the  center  of  A(x,x')  among  the 


Table  1. — HIRS  channel  spectral  characteristics 


Channel  Center 

C02  band  no<  frequency  (cm  x) 

15  ym  1  668.5 


680.0 


690.0 


703.0 


716.0 


733.0 


749.0 


4.3  ym  11  2190.0 


12  2210.0 


13  2240.0 


14  2270.0 


15  2360.0 
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Figure  1. — The  derivative  of  transmittance  with  respect   to   the  logarithm  of 
pressure  for  the  15-\xm  CO     absorption  band. 
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Figure  2. — The  derivative  of  transmittance  with  respect   to  the  logarithm  of 
pressure  for  the  4.3-\im  CO     absorption  band. 
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Figure   3. — Radiative  transfer  kernels  for  the  15-\xm  and  4.3~\\m  CO     absorption 

bands  vs.   pressure. 
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Figure  4. — Averaging  kernels  at  5  mb  for  the   4,3-um  band,   15-\xm  band  and 
4.3-\xm  +  15-\im  bands  plotted  as  functions  of  atmospheric  pressure. 
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Figure  5. — Averaging  kernels  at  10  mb  for  the  4.3-V-m  band,   15-]im  band,   and 
4.3-um  +  15-ym  bands  plotted  as  functions  of  atmospheric  pressure. 
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Figure  6. — Averaging  kernels  at  20  mb  for  the   4.3-um  band,   15-\im  band,   and 
4.3-]im  +  15-\xm  bands  plotted  as  functions  of  atmospheric  pressure. 
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Figure  7. — Averaging  kernels  at   50  mb  for  the  4.3-\xm  band,   15-\xm  band,   and 
4.3-\im  +  15-\im  bands  plotted  as  functions  of  atmospheric  pressure. 
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Figure  8. —Averaging  kernels  at  70  mb  for  the  4.3-\im  band,   15-\ym  band,   and 
4.3-\im  +  15-\im  bands  plotted  as  functions  of  atmospheric  pressure. 
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Figure  9. — Averaging  kernels  at  100  mb  for  the  4.3-]xm  band,   15-]sm  band,   and 
4.3-\xm  +  15-\im  bands  plotted  as  functions  of  atmospheric  pressure. 
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Figure  10. — Averaging  kernels  at  200  mb  for  the   4.3-um  band,   I5-um  band,   and 
4.3-u/n  +  15-\xm  bands  plotted  as  functions  of  atmospheric  pressure. 
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Figure  11. — Averaging  kernels  at   500  mb  for  the  4.3-\im  band,   15-\im  band,   and 
4.3-\xm  +  15-\im  bands  plotted  as  functions  of  atmospheric  pressure. 
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Figure  12. — Averaging  kernels  at   700  mb  for  the  4.3-\xm  band,   15-\im  band,   and 
4.3-uzn  +  15-\im  bands  plotted  as  functions  of  atmospheric  pressure. 
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three  cases.   Finally,  the  center  and  the  resolving  of  A(x,x')  for  the  above- 
mentioned  three  cases  are  computed  and  plotted  in  figures  13  and  14,  respec- 
tively. 

4.   DISCUSSION 

Qualitatively,  the  curves  of  figures  4  to  12  indicate  that  the  12  channels 
of  the  combined  4.3-um  and  15-um  CC>2  bands  has  better  vertical  resolution 
of  temperature  profiles  than  that  of  the  7  channels  of  the  15-um  band, 
uniformly  for  all  p,  which  in  turn  has  better  vertical  resolution  than  that 
of  the  5  channels  of  the  4.3-um  band.   This  is  substantiated  by  the  center 
vs.  p  curves  in  figure  13  and  the  resolving  length  w(p)  vs.  p  curves  in 
figure  14.   It  is  found  that  the  addition  of  the  5  channels  of  the  4.3-um 
band  to  the  15-um  band  improves  the  resolving  length  of  the  averaging  kernel 
for  the  15-utii  band  alone  about  1  to  15%  in  the  troposphere,  15  to  43%  in 
the  stratosphere  except  in  the  range  of  3.5  to  8.5  mb,  and  15%  at  the  tropo- 
pause.   Note  that  the  w(p)  vs.  p  curve  for  the  15-um  band  is  comparable  to 
that  of  Conrath  (1971).   Furthermore,  figure  13  shows  that  in  general  the  cen- 
ter of  the  averaging  kernel  for  the  combination  of  4.3-um  and  15-um  bands 
is  closer  to  the  level  to  which  the  kernel  pertained  than  the  other  two  cases. 
This  is  particularly  true  in  the  upper  stratosphere,  for  example  in  the  layer 
of  p  <  10  mb.   Although  the  resolving  length  of  the  15-um  band  is  shorter  than 
that  of  the  combined  4.3-um  and  15-um  bands  in  the  range  of  3.5  to  8.5  mb , 
the  large  deviation  of  the  center  of  the  averaging  kernel  for  the  15-um  band 
from  the  level  to  which  the  kernel  pertains  does  away  with  any  advantage 
of  shorter  resolving  length.   In  our  opinion,  these  facts  indicate  the 
uniform  superiority  of  the  combination  of  4.3-um  and  15-um  bands  over  the 
other  two  cases  for  all  p  in  the  vertical  resolution  of  temperature  profiles. 

Figure  14  indicates  that  the  resolving  length  for  the  4.3-um  C0~  channels 
is  poor  in  the  range  of  8  to  60  mb.   This  is  mainly  because  of  the  lack  of 
weight  concentration  of  K[x,v^]'s  in  that  particular  range  (figure  4). 
Unfortunately,  the  method  of  Backus  and  Gilbert  will  not  yield  accurate 
results  under  this  circumstance.   To  obtain  more  accurate  results,  one  has 
to  experiment  with  various  J-f unctions  other  than  that  of  eq  (12). 

It  is  known  that  the  instrument  noise  will  contaminate  the  above  results 
(Backus  and  Gilbert,  1970).   The  general  effect  of  the  noise  here  is  be- 
lieved to  be  the  same  as  that  in  the  analysis  of  Conrath  (1971). 
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